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Abstract. Theoretical approaches to explaining the observed acceleration of the 
universe are reviewed. We briefly discuss the evidence for cosmic acceleration, and 
the implications for standard General Relativity coupled to conventional sources of 
energy-momentum. We then address three broad methods of addressing an accelerating 
universe: the introduction of a cosmological constant, its problems and origins; the 
possibility of dark energy, and the associated challenges for fundamental physics; and 
the option that an infrared modification of general relativity may be responsible for 
the large-scale behavior of the universe. 
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1. Introduction 

The development of General Relativity almost 90 years ago provided not only a new 
way to understand gravity, but heralded the dawn of an entirely new branch of science 
- cosmology. Of course, a discussion of the visible contents of the night sky was not a 
new development - these objects and phenomena have forever fascinated civilizations. 
However, a coherent theoretical framework within which to explore their distribution 
and evolution was lacking before Einstein's crucial insight - that space and time 
are themselves players in the cosmic drama. Bubble's 1921 observation that distant 
galaxies are receding from ours at speeds proportional to their distances, in agreement 
with Friedmann's corresponding solution to General Relativity, then provided the first 
experimental confirmation of this new science. 

Almost a century has passed since the beginning of this era, and in the intervening 
years increasingly accurate predictions of this model of the cosmos, supplemented 
only by the presence of a dark matter component, have been confronted with, and 
spectacularly passed, a host of detailed tests - the existence of the Cosmic Microwave 
Background Radiation (CMB); the abundances of the light elements through Big Bang 
Nucleosynthesis (BBN); the formation of structure under gravitational instability; the 
small temperature anisotropics in the CMB; the structure of gravitational Icnsing maps; 
and many more. Many of these tests are highly nontrivial and provide remarkable 
support for the overall big bang model. 

All this held true until approximately a decade ago, when the first evidence 
for cosmic acceleration was reported by two groups using the lightcurves of type la 
supernovae to construct an accurate Hubble diagram out to redshift greater than 
unity. As we will describe in some detail, the natural expectation is that an expanding 
universe, evolving under the rules of GR, and populated by standard matter sources, 
will undergo deceleration - the expansion rate should slow down as cosmic time unfolds. 
Amazingly, that initial data, supplemented over the last decade by dozens of further, 
independent observations, showed that the universe is speeding up - what we observe is 
not deceleration, but cosmic acceleration! 

The goal of this article is to provide an overview of what is meant by cosmic 
acceleration, some proposed theoretical approaches to understanding the phenomenon, 
and how upcoming observational data may cast crucial light on these. We have not 
attempted to be completely comprehensive, but have tried to provide a flavor of the 
current research landscape, focusing on the better known approaches, and highlighting 
the challenges inherent in constructing a fundamental physics understanding of the 
accelerating universe. There are a number of excellent recent review articles on 
dark energy, which the reader might find as a useful complement to the present 
paper [1, 2, 3, 4]. 

The article is organized as follows. In the next section we provide a brief 
introduction to the relevant cosmological solutions to General Relativity, with perfect 
fiuid sources. In section 3 we review the observational evidence, primarily from 
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observations made during the last decade or so, for an accelerating universe. In 
section 4 we then discuss the minimal possibility that Einstein's cosmological constant 
is responsible for late-time cosmic acceleration, and discuss the challenges the required 
magnitude of this parameter raises for particle physics. In 5 we then explore a class of 
dynamical dark energy models, in which the universe accelerates due to the evolution of a 
new component of the cosmic energy budget, under the assumption that the cosmological 
constant is negligibly small. In a similar vein, in 6 we then consider the possibility 
that the universe contains only standard model and dark matter sources of energy and 
momentum, but accelerates because GR is modified at large distances, leading to new 
self-accelerating solutions. Finally, before concluding, in 7 we then consider the challenge 
of distinguishing among these possible explanations for cosmic acceleration, and discuss 
how upcoming observational missions may help us to better understand this phenomena. 

A note on conventions. Throughout this article we use units in which h = c = 1, 
adopt the (—,+,+,+) signature and define the reduced Planck mass by = — 
SttG. 



ds^ = -dt^ + a^{t) 



2. Essential Elements of Background Cosmology 

Our goal here is to lay out the bare minimum for understanding the background 
evolution of the universe. By background in this context, we mean the dynamics 
pertaining to the homogeneous and isotropic description of spacetime, valid on the 
very largest scales, without reference to spatial perturbations of either the metric or 
matter fields^. 

The most general homogeneous and isotropic metric ansatz is the Priedmann, 
Robertson- Walker (FRW) form 

_^ + r\de^ + sin'edcp^)] , (1) 

where k describes the curvature of the spatial sections (slices at constant cosmic time) 
and a{t) is referred to as the scale factor or the universe. Without loss of generality, we 
may normalize k so that A; = +1 corresponds to positively curved spatial sections (locally 
isometric to 3-sphcrcs); A; = corresponds to local flatness, and k = —1 corresponds 
to negatively curved (locally hyperbolic) spatial sections. These local definitions say 
nothing about the global topology of the spatial sections, which may be that of the 
covering spaces - a 3-sphere, an infinite plane or a 3-hyperboloid - but it need not be, 
as topological identifications under freely-acting subgroups of the isometry group of each 
manifold are allowed. As a specific example, the k — spatial geometry could apply 
just as well to a 3-torus as to an infinite plane. 

Note that we have not chosen a normalization such that oq = 1. We are not free 
to do this and to simultaneously normalize \k\ = 1, without including explicit factors 
of the current scale factor in the metric. In the flat case, where k = 0, we can safely 
choose Oo = 1. 

I We follow closely the brief discussion in [5] 
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At cosmic time t, the physical distance from the origin to an object at radial 
coordinate r is given by d{t) — a{t)r. The recessional velocity of such an object due to 
the expansion of the universe is then given by 

v{t) = H{t)d{t) , (2) 

where H{t) = d/a (an overdot denotes a derivative with respect to t) is the Hubble 
parameter, the present day value of which we refer to as the Hubble constant Hq. 

The Hubble parameter is a useful function through which to parametrize the 
expansion rate of the universe. As we shall encounter later, it is convenient to define 
a second function, describing the rate at which the expansion rate is slowing down or 
speeding up - the deceleration parameter, defined as 

m - -| . (3) 

2.1. Dynamics: The Friedmann Equations 

The FRW metric is merely an ansatz, arrived at by requiring homogeneity and isotropy 
of spatial sections. The unknown function a{t) is obtained by solving the differential 
equations obtained by substituting the FRW ansatz into the Einstein equation 

R^, - ^Rg^, = SttGT^, . (4) 

The energy-momentum tensors T^j, describes the matter content of the universe. It 
is often appropriate to adopt the perfect fluid form for this 

T^,.^ip + p)U^U,+pgf,u , (5) 

where is the fluid four- velocity, p is the energy density in the rest frame of the fluid 
and p is the pressure in that frame. 

Substituting (1) and (5) into (4), one obtains the Friedmann equation 

and 



Here i indexes all different possible types of energy in the universe. The Friedmann 
equation is a constraint equation, in the sense that we are not allowed to freely specify 
the time derivative d; it is determined in terms of the energy density and curvature. 
The second equation an evolution equation. 

We may combine (6) and (7) to obtain the acceleration equation 

---^Y.^Pi + ^Pi) , (8) 

i 

which will prove central to the subject of this review. 



Approaches to Understanding Cosmic Acceleration 



5 



In fact, if we know the magnitudes and evolutions of the different energy density 
components pi, the Friedmann equation (6) is sufficient to solve for the evolution 
uniquely. The acceleration equation is conceptually useful, but rarely invoked in 
calculations. 

The Friedmann equation relates the rate of increase of the scale factor, as encoded 
by the Hubble parameter, to the total energy density of all matter in the universe. We 
may use the Friedmann equation to define, at any given time, a critical energy density, 

for which the spatial sections must be precisely flat (A; = 0). We then define the density 
parameter 

^^total = - , (10) 

Pc 

which allows us to relate the total energy density in the universe to its local geometry 
via 

^^total ^l^k = (11) 

It is often convenient to define the fractions of the critical energy density in each different 
component by 

= - . (12) 

Pc 

Energy conservation is expressed in GR by the vanishing of the covariant divergence 
of the energy-momentum tensor, 

V^T'^^ = . (13) 

For the FRW metric (1) and a perfect-fluid energy-momentum tensor (5) this yields a 
single energy-conservation equation, 

p + 3H{p + p)^0 . (14) 

This equation is actually not independent of the Friedmann and acceleration equations, 
but is required for consistency. It implies that the expansion of the universe (as specified 
by H) can lead to local changes in the energy density. Note that there is no notion of 
conservation of "total energy," as energy can be interchanged between matter and the 
spacetime geometry. 

One final piece of information is required before wc can think about solving our 
cosmological equations: how the pressure and energy density arc related to each other. 
Within the fluid approximation used here, we may assume that the pressure is a single- 
valued function of the energy density p = p{p)- It is often convenient to define an 
equation of state parameter, w, by 

p = wp . (15) 
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This should be thought of as the instantaneous definition of the parameter w; it does not 
represent the full equation of state, which would be required to calculate the behavior 
of fluctuations. Nevertheless, many useful cosmological matter sources do obey this 
relation with a constant value of w. For example, w = corresponds to pressureless 
matter, or dust - any collection of massive non-relativistic particles would qualify. 
Similarly, w = 1/3 corresponds to a gas of radiation, whether it be actual photons 
or other highly relativistic species. 

A constant w leads to a great simplification in solving our equations. In particular, 
using (14), we see that the energy density evolves with the scale factor according to 

^(«) ^ a(t)3(l+^) • ^^^^ 

Note that the behaviors of dust {w — 0) and radiation {w — 1/3) are consistent with 
what we would have obtained by more heuristic reasoning. Consider a fixed comoving 
volume of the universe - i.e. a volume specified by fixed values of the coordinates, from 
which one may obtain the physical volume at a given time t by multiplying by a{tY. 
Given a fixed number of dust particles (of mass m) within this comoving volume, the 
energy density will then scale just as the physical volume, i.e. as a{t)~^, in agreement 
with (16), with w — Q. 

To make a similar argument for radiation, first note that the expansion of the 
universe (the increase of a{t) with time) results in a shift to longer wavelength A, or a 
redshift, of photons propagating in this background. A photon emitted with wavelength 
Ae at a time te, at which the scale factor is Og = a{te) is observed today {t — to, with 
scale factor ao = a{to)) at wavelength Ao, obeying 

+ (17) 

Ae fle 

The redshift z is often used in place of the scale factor. Because of the redshift, the 
energy density in a fixed number of photons in a fixed comoving volume drops with the 
physical volume (as for dust) and by an extra factor of the scale factor as the expansion 
of the universe stretches the wavelengths of light. Thus, the energy density of radiation 
will scale as a(i)~^, once again in agreement with (16), with w — 1/3. 

Thus far, we have not included a cosmological constant A in the gravitational 
equations. This is because it is equivalent to treat any cosmological constant as a 
component of the energy density in the universe. In fact, adding a cosmological constant 
A to Einstein's equation is equivalent to including an energy-momentum tensor of the 
form 

T,.^-:^g,.. (18) 

This is simply a perfect fluid with energy-momentum tensor (5) with 
A 

Pa = - Pa , (19) 
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so that the equation-of-state parameter is 

(20) 

This imphes that the energy density is constant, 

Pa = constant . (21) 

Thus, this energy is constant throughout spacetime; we say that the cosmological 
constant is equivalent to vacuum energy. 

Similarly, it is sometimes useful to think of any nonzero spatial curvature as yet 
another component of the cosmological energy budget, obeying 

3A; 

^^""^^ ~ ~ 8^ 

^'^"^^ = 8^ ' ^^^^ 

so that 

Wcurv = -1/3 . (23) 

It is not an energy density, of course; Pcurv is simply a convenient way to keep track of 
how much energy density is lacking, in comparison to a flat universe. 

2.2. Flat Universes 

Analytically, it is much easier to find exact solutions to cosmological equations of motion 
when A; = 0. We are able to make full use of this convenience since, as we shall touch on 
in the next section, the combined results of modern cosmological observations show the 
present day universe to be extremely spatially flat. This can be seen, for example, from 
precision observations of the cosmic microwave background radiation and independent 
measures of the Hubble expansion rate. 

In the case of flat spatial sections and a constant equation of state parameter w, 
we may exactly solve the Friedmann equation (16) to obtain 

a{t) = ao i-j , (24) 

where oq is the scale factor today, unless ■u; = — 1, in which case one obtains a{t) oc e^*. 
Applying this result to some of our favorite energy density sources yields table 1. 

Note that the matter- and radiation-dominated flat universes begin with a = 0; 
this is a singularity, known as the Big Bang. We can easily calculate the age of such a 
universe: 

Unless w is close to —1, it is often useful to approximate this answer by 

to ~ H^' ■ (26) 

It is for this reason that the quantity Hq^ is known as the Hubble time, and provides a 
useful estimate of the time scale for which the universe has been around. 
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Type of Energy 



p(a) 



a{t) 



Dust 

Radiation 

Cosmological Constant 



constant 



a 



4 




Table 1. A summary of the behaviors of the most important sources of energy density 
in cosmology. The behavior of the scale factor applies to the case of a flat universe; 
the behavior of the energy densities is perfectly general. 



3. Observational Evidence for Cosmic Acceleration 

An important breakthrough in cosmology occurred in the late '90s with the 
measurements of type la Supernovae (SNela) by the High-Z Supernova team [6] and 
the Supernova Cosmology Project [7]. As we review in this section, the data from 
these surveys, as well as from complementary probes, provide strong evidence that the 
universe has recently entered a phase of accelerated expansion. 

3.1. Type la Supernovae 

Type la supernovae are stellar explosions that occur as a white dwarf, onto which mass 
is gradually accreting from a companion star, approaches the Chandrasekhar limit [8] 
(of about 1.4 solar masses) and the density and temperature in its core reach the ignition 
point for carbon and oxygen. This begins a nuclear flame that fuses much of the star up 
to iron. They are extremely bright events, with luminosities a significant fraction of that 
of their host galaxies during the peak of their explosions. Therefore they arc relatively 
easy to detect at high redshift (2; ~ 1). This specific type of supernova is characterized 
by the absence of a hydrogen line (or, indeed, a helium one) in their spectrum (which 
is typical of type I SNe) and instead by the presence of a singly-ionized silicon line at 
615nm, near peak light. 

Despite a significant scatter, type la supernovae peak luminosities have been found 
to be very closely correlated with observed differences in the shapes of the light curves: 
dimmer SNe are found to decline more rapidly after maximum brightness, while brighter 
SNe decline more slowly [9, 10, 11]. This difference seems to be traceable to the amount 
of ^^Ni produced in the explosion; more nickel implies both a higher peak luminosity 
and a higher temperature, and thus opacity, leading to a slower decline. 

After one adjusts for the difference in the light-curves, the scatter in peak luminosity 
can be reduced to 15%. In this sense, SNela are referred to as "standardizable 
candles" and are very good candidates for distance indicators, since one would expect 
any remaining difference in their peak luminosity to be due to a difference in distance. 
Another important aspect of this uniformity is that it provides standard spectral and 
light curve templates that offer the possibility of singling out those SNe that deviate 
slightly from the norm [12]. 

After the necessary adjustments, all SNela have the same absolute magnitude M, 
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and any difference in tfieir apparent magnitude m is attributable to a difference in tlieir 
distance. Tfie apparent and intrinsic magnitudes are related via the luminosity distance 



where T^' is a correction for the shifting of the spectrum [6, 7, 9, 11, 13, 14] necessary 
because only a part of the spectrum emitted is actually observed. In an expanding 
universe, the expression of the luminosity distance as a function of redshift is 



By independently measuring dL{z) and z, one can then constrain the expansion history 
H{z). Following pioneering work reported in [15], the High-Z Supernova team [6] and 
the Supernova Cosmology Project [7] measured the apparent magnitudes of many SNela 
(of redshift z < 1) and directly determined their distances. They then compared these 
distances to those inferred from the redshifts of the host galaxies (measured from the 
spectra of the galaxies when possible, otherwise the spectra of the SNe themselves). The 
most distant SNe appeared dimmer than expected in a universe currently dominated by 
matter, and a careful comparison with low redshift supernovae allowed one to rule out 
the possibility that this dimming is due to intervening dust. Assuming a fiat universe 
described by GR, and homogeneous and isotropic on large scales, the data are best fit by 
a universe which has recently entered a phase of accelerated expansion, i.e. a universe 
for which the deceleration parameter (3) is currently negative {qo ~ — 1). The simplest 
model fitting the data is a universe in which matter accounts for only about a quarter of 
the critical density, while the remaining 70% of the energy density is in the cosmological 
constant A. This model is commonly referred to as the ACDM model. In fact, the SNe 
data by themselves allow a range of possible values for the matter and cosmological 
constant density parameters {Qy^ and respectively). However, as we will discuss in 
the following subsections, we can use what we know from other observations, such as the 
Cosmic Microwave Background (CMB) and the large scale structure (LSS), to further 
constrain these parameters. 

In more recent years, distant SNe with redshifts up to z ~ 1.7 have been observed 
by the Hubble Space Telescope (HST) [17, 18]. These observations show that the trend 
toward fainter SNe at moderate redshifts has reversed; therefore they play a key role 
in disregarding dust as a plausible explanation of the dimming of intermediate SNe, 
and in firmly establishing [19, 20] acceleration. More generally, SNe at lower and 
higher redshift are both important in the study of the expansion history. Nearby SNe 
{z < 0.3) are useful for defining characteristics of Type la SNe, understanding the 
explosion and exploring systematics; all combined, they are important for establishing 
distance indicators and are currently [21, 22] being measured, with more to come in 
the future [23, 24]. Intermediate redshift SNe (0.3 < < 0.8) measure the strength 
of cosmic acceleration and are being measured by the ESSENCE project [25], the 



dL 




(27) 




(28) 
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Figure 1. Hubble diagram for Type la Supernovae, plotting the effective magnitude 
TTiB versus redshift z (Knop et al. [16]). The solid line is the best-fit cosmology and 
corresponds to a universe with il^ = 0.25 and fi^ = 0.75. 

SuperNova Legacy Survey (SNLS) [26] and the Sloan Digital Sky Survey-II (SDSS- 
II) Supernova Survey [27]. Finally, distant SNe {z > 0.8) are important to break the 
degeneracy among the cosmological parameters and they are currently being measured 
by the Higher-z Supernova Search Team (HZT) [28] and the GOODS team [29]. 

In the near future, the space-based Joint Dark Energy Mission (JDEM) [30], 
(a wide-field optical-infrared telescope), will offer precision measurements with the 
aim of determining the nature of cosmic acceleration. This mission may incorporate 
the Baryon Acoustic Oscillation (BAO), Supernovae (SNe) and Weak Lensing (WL) 
techniques. Also, upcoming and future surveys such as the earth-based Dark Energy 
Survey (DES) [31] and the Large Synoptic Survey Telescope (LSST) [32] will provide 
observations useful for constraining the expansion history of the universe. 

3.2. Complementary Probes 

Supernovae are the first and most direct evidence for the late-time acceleration of the 
universe. However, independent evidence comes from complementary probes. Data 
from different observations are important not only to confirm the results from SNe, 
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but also to help break some of the degeneracy in the cosmological parameters, such as 
the density parameters and Q^. As we have mentioned, the measurement of the 
expansion history via SNe is an example of the use of standard candles to infer the 
luminosity distance as a function of the redshift d^lz). Alternatively, one may employ 
standard rulers, extracting the angular diameter distance dA{z) (given by the ratio of 
the object's actual size to its angular size) and then H{z). (Both methods are commonly 
referred to as geometrical probes). The expansion history can further be constrained 
with the growth of structure, as we will see shortly. Finally, some methods are based 
on inferring the age of the universe from stellar ages. In what follows we review several 
complementary probes of the expansion history of the universe. 

The CMB is the radiation that reaches us from the surface of last scattering - that 
is, from the time at which photons decoupled from the thermal bath and started to 
free stream to us. It is in the form of an almost isotropic blackbody spectrum, with a 
temperature of approximately 2.7 K. The observed blackbody distribution, typical of 
sources that are in thermal equilibrium, is strong evidence for the theory of the Big 
Bang. At early times, when the universe was radiation dominated, photons were easily 
energetic enough to ionize hydrogen atoms. Therefore, the universe was filled with a 
charged plasma, the photon-baryon plasma, and was opaque. As the universe expanded, 
the temperature decreased and the photons redshifted, losing energy until they were 
no longer able to ionize hydrogen. At T ~ 0(1) cV, corresponding to a redshift of 
z ~ 1100, neutral hydrogen formed and the photons decoupled from the plasma. At 
early times, the number and energy densities of matter were high enough that the 
plasma constituents were in thermal equilibrium, ensuring a blackbody spectrum for 
the radiation distribution. This spectrum is maintained by photons after decoupling, 
just at lower and lower temperatures. Indeed, the effect of homogeneous and isotropic 
expansion on a blackbody distribution is that of preserving the blackbody structure 
with temperatures that redshift with the wavelength, T oc 1/a. 

As measured by the Cosmic Background Explorer (COBE) [33], the temperature 
distribution of the photons is not completely isotropic, with fluctuations of one part 
in 10^. These fluctuations contain a great deal of physical information about our 
universe, and the ongoing measurements from the Wilkinson Microwave Anisotropy 
Probe (WMAP) [34, 35, 36] have reached unprecedented precision in using it to constrain 
cosmological parameters. It is worth exploring in a little more detail how this is achieved. 

Matter in the universe today is highly inhomogeneous and clustered. Our 
understanding of how this large scale structure developed is that initially small density 
perturbations, in an otherwise homogeneous universe, grew through gravitational 
instability into the objects we observe today. The basic observable in CMB physics 
is the temperature fluctuation on the sky, for which we measure the 2-point correlation 
function. The Fourier counterpart of the correlation function, the power spectrum, is 
what is commonly shown in plots. As can be seen in Fig. 2 , there are two major 
components of a typical CMB power spectrum: the Sachs- Wolfe plateau at large scales 
and a series of peaks and troughs at smaller scales. These both carry important 
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Figure 2. The Cosmic Microwave Background angular power spectrum, representing 
the two-point correlation function of temperature fluctuations (in momentum space) 
on the sky as a function of angular separation, [courtesy of NASA/WMAP Science 
team]. The solid line is the best fit ACDM model. On larger scales, i.e. at lower 
multipoles I, we notice the Integrated Sachs- Wolfe plateau, at I ~ 220 we see the first 
acoustic peak and on smaller angular scales the subsequent acoustic peaks 



information about the source of cosmic acceleration. The fluctuations we observe 
on the largest scales are essentially the primordial ones, since the largest scales have 
only recently entered the horizon and therefore little causal physics has affected them. 
On smaller scales we instead observe fluctuations which entered the horizon before 
decoupling and have undergone acoustic oscillations characteristic of the photon-baryon 
fluid. 

Let us focus on the first peak of these oscillations. The acoustic oscillations are 
damped oscillations of the photon-baryon plasma under the competing effects of gravity 
and radiation pressure. The maximum amplitude of the oscillations is expected on 
fluctuations that entered the horizon just before decoupling. To see this, consider an 
overdense region of size R, which therefore contracts under self-gravity over a timescale 
R (since we have set c = 1). If R ^ Hj^^ then the region will not have had time to 
collapse over the lifetime of the universe at last scattering. If R <^ Hj^^ then collapse will 
be well underway at last scattering, matter will have had time to fall into the resulting 
potential well and result in a rise in temperature. This, in turn, gives rise to a restoring 
force from photons which acts to damp out the inhomogeneity. Clearly, therefore, the 
maximum anisotropy will be on a scale which has had just enough time to equilibrate, 
i.e. R ~ Hf^^. This means that we expect to see a peak in the CMB power spectrum 
at an angular size corresponding to the horizon size at last scattering. Since we know 
the physical size of the horizon at last scattering, this provides us with a ruler on the 
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Figure 3. Observational constraints in the 17^ — 57° plane [37]. The contours represent 
the 68.3%, 95.4% and 99.7% confidence level from CMB, BAO and SNela and their 
combination, with different colors corresponding to different data sets. The vertical 
green contours, centered around ~ 0.3, represent BAO constraints from the SDSS 
galaxy survey [38] ; the orange narrow diagonal contours correspond to constraints from 
WMAP-5 year observations of the CMB anisotropics (including a prior on the Hubble 
parameter) [39]. Finally the wide blue diagonal contours represent constraints from 
supernovae Union Set [37] and the gray-scale contours the combined constraints from 
BAO, SNela and CMB. 



sky. The corresponding angular size will then depend on the spatial geometry of the 
universe. For a flat universe {k = 0, fl^^^ = 1), we expect a peak at / ~ 220. As can 
be seen in Fig. 2, the first peak indeed appears at this angular scale, providing strong 
evidence for the spatial flatness of the universe. As represented in Fig. 3, the CMB 
strongly constrains the sum of the matter and dark energy densities, offering a probe 
which is complementary to the SNe data in the — space. 

However, this is just one piece of the valuable information that can be extracted 
from the CMB. In general, the CMB power spectrum provides high precision constraints 
on essentially all of the cosmological parameters [36]. The overall CMB spectrum, as 
well as the matter power spectrum as measured from LSS, are best fit by models which 
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Figure 4. Schematic representation of the Integrated Sachs- Wolfe effect (ISW) on 
CMB photons. On the left side of the figure, a photon travels through a constant 
potential, gaining and losing the same amount of energy while passing through the 
potential well. In this case, the outgoing wavelength corresponds to the incoming one, 
and there is no net gain (or loss) of energy. On the right side of the figure, a photon 
travels through a shallowing potential. As the potential well from which the photon 
exits is shallower than the one it entered, there is an overall gain in energy and the 
outgoing wavelength is shorter than the incoming one. 



undergo late-time accelerated expansion. 

Information from the CMB concerning the total energy density of the universe can 
be combined with data about the distribution of matter, from which one extracts In 
addition, the CMB power spectrum itself offers a way to constrain the matter density via 
the spacing of the peaks of the acoustic oscillations. Complementary strong constraints 
on the matter density come from the observed distribution of matter in the universe. 
In particular, the turnover scale in the matter power spectrum is very sensitive to 
and the abundance of clusters of galaxies constrain the as parameter, which encodes 
information about and the fractional density perturbation. 

The WMAP data together with LSS data from the Sloan Digital Sky Survey 
(SDSS) [40] yield the value fi^ = 0.25 ± 0.03 (la, in the context of a flat ACDM 
model) [41]. Combining these with the supernovae data, we can pin down a narrow 
region m space corresponding to fi^ ~ 1 — ~ 0.27, as shown in Fig. 3. 

Gravitational lensing - the distortion of light paths by the inhomogeneous metric 
around matter in the universe - is also an important complementary probe as it offers a 
direct handle on the mass distribution, without relying on bias factors to convert from 
visible to total matter. Of particular importance in this context is the weak gravitational 
lensing of distant galaxies by intervening foreground mass overdensities, dubbed cosmic 
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Figure 5. The large scale redshift-space correlation function of the SDSS LRG 
(luminous red galaxies) sample, from [38]. The magenta line shows a pure CDM 
model (ri^/i^ = 0.105) which lacks the acoustic peak. The green, red and blue lines 
represent models with a baryon fractional density of fi^/i^ = 0.024 and a CDM density 
of, respectively, fi^^^ = 0.12,0.13,0.14. The bump at 100 h~^Mpc corresponds to the 
acoustic peak. Details about the data set and error bars can be found in [38]. 

shear. This offers a direct measurement of mass ffuctuations in the universe and their 
redshift evolution. The statistic associated with cosmic shear, the shear angular power 
spectrum, is sensitive to dark energy via the expansion history as well as via the growth 
of structure. From precise measurements of the cosmic shear and its cross-correlation 
with galaxies it is possible to extract geometrical probes of dark energy [42, 43] as well 
as an estimate of the matter density parameter [44]. In section 7 we will discuss in more 
detail the constraints on the growth of perturbations from weak lensing tomographic 
surveys. 

Further useful constraints on the matter density can be derived from measurements 
of the fraction /gas of X-ray emitting gas to total mass in galaxy clusters. This ratio 
is a good indicator of the overall baryon fraction in the universe, and depends on the 
distance to the cluster, offering a probe of the expansion history. The ratio of the gas 
(baryons) to total mass is found to be about 11%, which yields a cosmic baryon fraction 
of about 15% [45, 46], after estimated corrections. Combining these measurements with 
the baryon density inferred very precisely from Big Bang Nucleosynthesis, one can thus 
derive an estimate of the matter density parameter. 

The large scale plateau of the CMB power spectrum offers a useful probe of the 
expansion history. In a matter dominated universe, when the largest fluctuations enter 
the horizon, the gravitational potentials are constant, and therefore do not imprint any 
additional anisotropy on the photons that traverse them. Indeed, a photon entering a 
potential well gains energy, while a photon climbing out of a potential loses energy. 
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as represented in Fig. 4; if the potential is constant in time, there is no overall 
change in the energy of the photon. However, if matter is not dominant and the 
late universe is undergoing a phase of accelerated expansion, then the gravitational 
potentials are expected to evolve. As a consequence, photons traveling through a 
shallowing (deepening) potential well will be blueshifted (redshifted) , as they gain more 
energy entering the well than the energy they lose climbing out. This is known as the 
Integrated Sachs- Wolfe effect (ISW). However, this is a small contribution to the overall 
CMB power spectrum on the largest scales, approximately 10% of the total signal, and 
therefore is very hard to extract. One proceeds by cross-correlating the CMB power 
spectrum with the galaxy one, as the ISW signal is expected to be strongly correlated 
with the distribution of matter around us. At present, the ISW effect has been measured 
out to redshift 2; ~ 1.5 [47, 48] at a 4(T significance level, and constitutes a promising 
probe of dark energy that is likely to play an important role in distinguishing among 
candidate models of cosmic acceleration, as we will discuss in section 6. However, it 
must be noted that the ISW effect will always be hmited by a low signal to noise ratio. 

The acoustic oscillations that we observe in the spectrum of the CMB also leave 
a characteristic imprint on the matter power spectrum [49, 50, 51, 52], commonly 
referred to as Baryon Acoustic Oscillations (BAO). These imprints provide a new 
opportunity to measure distances at different redshifts. Indeed, the length scale that was 
imprinted by the sound waves, referred to as the acoustic scale, persists to the present 
time in the clustering of matter at a characteristic comoving length of 100 h^^ Mpc. 
The latest generation of galaxy redshift surveys is able to probe the length scales 
required to make a precision measurement of the BAO signal, and this method has 
therefore been proposed as a new geometrical probe to constrain the physics of cosmic 
acceleration [53, 54, 55, 56, 57]. The acoustic scale corresponds to a preferred length 
scale on the surface of last scattering, i.e. the distance that a sound wave propagating 
from a point source at the end of inflation would have traveled before decoupling. In a 
redshift survey, this length scale can be measured both along and across the line of sight, 
allowing the observer to reconstruct the Hubble parameter and the angular diameter 
distance as a function of redshift, respectively H{z) and (/^(-z). The BAO have been 
detected in the SDSS data [38] data at a redshift of approximately z ~ 0.35 and in 
a combination of the Two Degree Field Galaxy Redshift Survey (2dFGRS) and SDSS 
data [58]. 

From the recent WMAP 5-year results [36], the combined constraints from BAO, 
CMB and SNela on the density parameters are: = 0.0456 ± 0.0015, Q^/i^ ^ 
0.13581^;^^^^ and Ql = 0.725 ± 0.015 to 68% C.L.. 

4. The Cosmological Constant 

The natural tendency of a universe filled with regular matter is to decelerate, as the 
mutual gravitational attraction fights against any initial expansion. This was clear even 
to Einstein, for whom this fact was an obstacle to obtaining a static universe, rather 
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than an accelerating one. Indeed, as Einstein realized almost immediately, there is an 
extension of the field equations of general relativity, consistent with general covariance 
and second order field equations, that can battle contraction and ultimately lead to 
acceleration - the cosmological constant. 

The cosmological constant A is the coefficient of the term in the field equations 
that is proportional to the metric. It appears in the action for gravity as 



where T^i, is the energy-momentum tensor for mass-energy sources in the universe. 

As we have mentioned earlier, for cosmological purposes, we may treat the 
cosmological constant as a component of the energy density in the universe, since adding 
a cosmological constant to the Einstein equation is equivalent to including an energy- 
momentum tensor of a perfect fluid with equation-of-state parameter is uja = ~1- 

As we have also seen, the cosmological constant is a prime candidate for the agent 
causing cosmic acceleration and, subject to ever more accurate data, provides a rather 
good fit to all known cosmological observations. Furthermore, in classical general 
relativity the cosmological constant is merely a constant of nature, with dimensions 
of [length]"^, in the same sense that the Planck mass in nothing but a dimensionful 
constant of the theory. In this setting, it is meaningless to enquire about the value 
of the constant, rather it is just something that we determine through experiment. 
However, when we begin to think about the quantum nature of matter it becomes clear 
that the cosmological constant poses a deep puzzle for fundamental physics. 

Particle physics, described by quantum field theory, is a remarkably accurate 
description of nature, and one of the best tested theories in all of science. Within 
this theory, we can attempt to calculate the expected contribution to the cosmological 
constant from quantum fluctuations in the vacuum. Taking the Fourier transform of a 
free (noninteracting) quantum fleld, each mode of flxed wavelength is a simple harmonic 
oscillator, which we know has zero-point (ground state) energy £'0 = jhw, where cu is 
the frequency of the oscillator. The net energy of these modes is then given by a 
divergent integral. Thus, our most naive expectation of the vacuum energy calculated 
within field theory is that it should be infinite! Perhaps, however, we have been too 
hasty in performing this calculation. After all, we have included modes of arbitrarily 
small wavelength, and this presumes that we can trust our descriptions of matter and 
gravity down to such tiny scales. Without a well-understood and calculable quantum 
theory of gravity, this is surely not reasonable. A better estimate of the vacuum energy 
can therefore be obtained by introducing a cutoff energy, above which we ignore the 
contributions of any modes, assuming that this will be justified by a more complete 
theory. This approach corresponds to assuming that effective fleld theory correctly 




(29) 



and in the resulting Einstein fleld equations as 




(30) 
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describe the appropriate limit of quantum gravity. In that sense, if a parameter has 
dimensions [mass]**, then we expect the corresponding mass parameter to be driven up 
to the scale at which the effective description breaks down. Thus, if classical general 
relativity is valid up to the Planck scale, we expect the vacuum energy to be given by 

^(theory) ^ ^4 _ (3^) 

Comparing this value to the value piac^^ required to explain cosmic acceleration, we 
obtain 

^(obs) ^ io-120^C«jeory) ^ (33) 

which is embarrassing, to say the least. It is perhaps more sensible to express the 
vacuum energy in terms of a mass scale as pvac = ^vao terms of which the value 
required to explain our observations is M^ac^'* ~ 10~^ eV, satisfying 

M^°^^) ~ 10-^°M^*^^°'^y) . (33) 

Nevertheless, this discrepancy of 30 orders of magnitude in energy is overwhelming, and 
is what is meant by the cosmological constant problem. 

One may add to this problem the following puzzling observation. The ratio of the 
vacuum and matter densities changes as the universe expands according to 

^ = -oca^ (34) 

Thus, only during a brief epoch of cosmic history is it possible for observers to witness 
the transition from matter domination to A domination, during which and flu are 
of the same order of magnitude. This is known as the coincidence problem. 

The issue of rehably calculating the cosmological constant, and finding a framework 
in which that calculation leads to a result dramatically different to the expected one 
has proven remarkably resistant to theoretical attack. It is fair to say that there are not 
currently any especially promising approaches. Nevertheless, there are two theoretical 
lines of research that are worth mentioning in this context. 

The first is supersymmetry (SUSY). Supersymmetry is a spacetime symmetry that 
relates particles of different spins (there exists a matching fermionic degree of freedom 
for every bosonic one, and vice-versa.) Since bosons and fermions make opposite sign 
contributions to the vacuum energy, in a truly supersymmetric theory the net vacuum 
energy should be zero, with cancellations occurring between, for example, the spin-1/2 
electron and its spin-0 partner, the "selectron" (with the same mass and charge as the 
electron.) 

It is, however, entirely obvious from experiment and observation, that the world 
is not supersymmetric - we would certainly know if there existed a scalar particle with 
the same mass and charge as an electron. Thus, if supersymmetry is realized in nature, 
it must be broken at some scale Afgusy which current bounds force to be greater than 
around a TeV. An analogous calculation of the vacuum energy in such a theory then 
yields 

Mvac ~ M,usy , (35) 
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with pvac = -^vac- That SUSY is a beautiful idea provides us with no clue as to the 
possible value of Mgusy However, independent of the cosmological constant problem, or 
any cosmological data, SUSY has been suggested as a solution to one of the outstanding 
problems of the standard model of particle physics - the hierarchy problem: why is the 
scale of electroweak symmetry breaking so much smaller than the scale of gravity. SUSY 
can help with this problem if the SUSY breaking scale is very close to current bounds, 
i.e. 

Msusy ~ 10^ GcV , (36) 

leading to the hope that we may find evidence for SUSY at the upcoming Large Hadron 
Collider (LHC). However, even this lower bound on the SUSY breaking scale leaves us 
with 

Mi:^'^ ~ lO-''M,^y , (37) 

which is still nowhere near acceptable. 

Nevertheless, should we find evidence for supersymmetry through particle physics 
experiments, the question of the strangely small value of the cosmological constant will 
become one of the SUSY breaking sector of the theory, and we will have made some 
progress, although how much is unclear. 

A second theoretical approach to the cosmological constant and cosmic acceleration 
problems is to consider the logical possibility that the vacuum energy is a feature 
of our local environment - the anthropic principle. This is an idea that may make 
logical sense if there exist many different possible values of physical parameters (and 
the cosmological constant in particular), and a mechanism in the early universe for 
sampling these possibilities, either over and over again in time, or in different regions of 
space. 

Clearly, only a certain range of values of the cosmological constant is consistent 
with the development of life like us (and there are a number of assumptions here), since 
if |A| is significantly larger than that required to yield the observed cosmic acceleration, 
then the universe would have recoUapsed (if Pvac were negative) , or expanded so quickly 
that galaxies couldn't form (if pvac were positive). Given this, some authors have 
argued [59, 60, 61, 62, 63, 64] that, under an assumption about the distribution of 
possible values of pvac, we should expect a value of A close to that we seem to have 
observed. 

The recent resurgence of interest in this possibility stems from the realization 
that within string theory there seems to exist a large number (more than 10^°*^) of 
possible vacuum states [65, 66, 67, 68, 69, 70]. These vacua arise from different 
compactification possibilities, and the roles played in them by various brane and 
gauge field configurations. The low energy effective theories constructed around these 
metastable vacua should have different local values of the vacuum energy (and other 
physical parameters.) An approach to addressing the cosmological constant problem 
then begins by asserting that it is highly likely that there must exist (and, eventually, 
hopefully identifying) one or more vacua with a cosmological constant (and perhaps 
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other parameters) that take the values one might think this reasonable given the number 

quoted above. 

The second step in constructing an anthropic argument around the string theory 
landscape involves arguing that there is a mechanism to sample, or populate, the whole 
distribution of possible vacua, in order to be able to claim that if a phenomenologically 
viable vacuum exists then some part of the universe will exist in that state. The 
mechanism of choice in the string theory landscape is eternal inflation. Inflation itself 
works to expand a small region of space to a size larger than the observable universe and, 
in general, models of inflation tend to be eternal, meaning that the universe continues 
to inflate in some regions even after inflation has ended in others [59, 71] . In this model 
many different domains may be inflated separately, and then small parts of the universe 
may cease to inflate and drop out into different final vacuum states. Thus, if there 
exist possible vacua with a cosmological constant of the right size to explain cosmic 
acceleration, parts of the universe that populate those vacua will cease inffating, reheat 
the universe, undergo standard cosmic evolution, and ultimately begin accelerating again 
at late times. 

This is a brief description of how the anthropic principle may be realized in the 
context of the string theory landscape. It is however, too early to know whether it can 
be realistically implemented, and even harder to say whether one will ever be able to 
falsify, and hence make science out of, such an approach. 

Nevertheless, this is the best developed argument regarding the possible small value 
of A, and perhaps for the value needed to explain cosmic acceleration. 

5. Dark Energy 

As we discussed in section 3, while the data clearly points to an accelerating universe, 
it does not single out a cosmological constant as the clear source of this acceleration. 
Let us assume, therefore, that the cosmological constant is zero, or at least so small as 
to be unable to account for the current period of accelerated expansion, and explore 
what other sources of mass-energy in the universe may be responsible (sec [72] for a 
discussion of logical possibilities, and [1] for a detailed review). Cosmologists have 
asked this question before, when confronted with the inflationary paradigm for solving 
the problems of early universe cosmology, and so we do not come the issue completely 
cold. 

The most faithful servant of the theoretical cosmologist is the trusty scalar field. 
It has served us dutifully not only through inflation, but also baryogenesis models, 
topological defect formation and some dark matter models, and so it is only natural 
that we turn to it now, in our epoch of need. Consider a real scalar field (p, with 
potential V{(j)), with Lagrangian density 

>C = ^^'^'^9^9,0 - y(0) , (38) 
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for which the energy-momentum tensor is given by 



]^<^Pd„(pd,(P + V{(P) 



(39) 



If the scalar field is homogeneously distributed in space, as we shall assume, then (39) 
takes the perfect fluid form, with energy density and pressure given by 
1 



1 



2r + n0), 



p = - V{<t^) , 

and we may then associate with it an instantaneous equation of state parameter 

„,=e^^. (40) 

As we learned when studying inflation, in the slow-roll regime, in which the scalar 
field equation of motion 

+ 3i/(/. + l^,^ = , (41) 

(where V^^ = dV/d(j)) becomes ~ ~3^^</" kinetic energy contributions to 
the pressure and the energy density become subdominant, and the equation of state 
parameter can become sufficiently negative to drive cosmic acceleration. 

Unlike inflation however, we need not require this period of acceleration to end, 
or the universe to be reheated after it. Rather, our requirements are that the energy 
scale associated with the potential be such that acceleration is happening at the current 
epoch, and that the energy density in the scalar field comprise the measured value of 
r^DE- Two of the most studied examples of this are inverse power law potentials 

and exponential potentials 

y(0)=M^exp(^-aA^ , (43) 

where M is some mass scale and a is a dimensionless parameter. 

In the latter case, it is well known that in the limit that the scalar field becomes so 
dominant that one can neglect contributions from other sources (such as matter), the 
scalar field and Priedmann equations can be solved exactly, to yield 

, 2, / M^a^ t 
(f){t) = - In 



a{t) = aot'/°' . 

Therefore, acceleration (a > 0) is achieved if < 2. 

Including matter changes this conclusion of course. One effect is that the actual 
functional form of the accelerating solution is altered from the simple power law one. 
However, a more interesting modification is the possible existence of tracker solutions 
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(in which the energy density in the scalar field approaches a fixed fraction of that in the 
other components, rather insensitively to the initial conditions), which provide a way 
to ameliorate the coincidence problem of cosmic acceleration. 

While these simple models provide a basic framework in which to obtain cosmic 
acceleration from a dynamical component, they are not without problems. Perhaps the 
largest issue is the extreme fine-tuning imposed for the same reasons as we described 
when discussing the cosmological constant. It is one thing to write down a potential, 
but quite another to provide a reason for the energy scales involved, and an argument 
why they should remain stable in any sensible particle physics implementation of the 
idea. 

To understand the fine-tuning, consider the mass of the quintessence field by 
expanding the potential about a point as V{(j)) = (l/2)m^0^. As wc discussed in the 
previous section, to match observations, the energy density in whatever component is 
responsible for cosmic acceleration must be ~ (10~^^ GeV)''. Assuming slow-roll, and 
that (0) ~ Mp, we may then equate these two to obtain 



Obtaining this kind of value from a particle physics model is an extreme challenge, 
consisting of two parts. First, it is not a simple matter to construct a model in which 
such low energy scales arise in the first place. Second, we face the issue confronting all 
scalar fields - how are we to protect an unnaturally small choice of mass scale against 
quadratic divergences when we compute quantum corrections to the couplings? 

As with inflation, one possibility is to choose a quintessence model in which a small 
mass is protected as a small parameter describing deviations from a point of enhanced 
symmetry in model space. This can be implemented [73] by making the quintessence field 
a Pseudo-Nambu-Goldstone Boson (PNGB) of a spontaneously broken shift symmetry. 
In the absence of soft breaking, the potential for such a field is a precisely fiat function 
of an angular coordinate. Thus, if the potential is lifted by a small mass term the exact 
symmetry keeping the potential fiat is only softly broken, and the potential takes the 
form 



where M is the overall mass scale, and / ~ Mp measures the range of variation of (f). 
Expanding this potential around = yields a field with a small mass (for appropriately 
small M), which is protected from quadratic divergences since the exact U{1) axion 
symmetry is restored in the f ^ oo limit. 

We shall not dwell on the simplest quintessence models in this review, since space 
is limited and our intention is to provide some idea as to more recent work in the field. 
The model wc have presented above admits a number of generalizations which have 
their own interesting features and signatures. 



~ lO"^'^ eV . 



(44) 




(45) 
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5.1. Phantom Models 

As we have described earlier, observations allow for the possibility that the equation of 
state parameter for dark energy obeys w < —1. However, thus far the models we have 
constructed do not venture into this parameter region. There is a good reason for this: 
in general relativity, it is conventional to restrict the possible energy-momentum tensors 
by imposing energy conditions in order to prevent instabilities of the vacuum or the 
propagation of energy outside the light cone. Standard examples, where in parentheses 
we have evaluated the condition for a perfect fluid source, are 



• The Weak Energy Condition: T^JH" > V timelike t" (p > and p + p > 0.) 

• The Null Energy Condition: T^^Ji^r > V null l^" {p + p> 0.) 

• The Dominant Energy Condition: T^U'^t" > V timehke ti", and Ti^^T^xtn^ < 



• The Null Dominant Energy Condition: T^Ji't > V null and T^''^^ is a 
non-spacelikc vector (as for the dominant energy condition, except that negative 
densities arc allowed if p = —p.) 

• The Strong Energy Condition: T^j^tH^ > ^T^xf^t^ V timelike t^ {p + p > and 
p + 3p> 0.) 

Since we are interested in an energy condition that implies stability, while allowing 
for the possibility of a cosmological constant, Garnavich [74] has advocated the use of 
the null dominant energy condition, boiling down to demanding that dynamical fields 
obey the dominant energy condition, while making an exception for the cosmological 
constant. Nevertheless, given what is allowed by the data, we may ask what happens 
if we allow for dark energy that violates the null dominant energy condition; i.e. obeys 
w < -1. 

The possibility of such phantom components has been suggested by many 
authors [75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85], and the simplest example [75] is 
provided by a scalar field x with negative kinetic and gradient energy 



obeying w < —1. 

While we will not describe the behavior of phantom fields in any detail, it is worth 
pointing out that, classically, the cosmology of such a system is relatively simple. For 
example, if the universe contains only dust and phantom matter. Then, if matter 
domination ceases at cosmological time t^, then the solution for the scale factor is 



(46) 



The associated equation of state parameter is 



Px_ v^(x) + |[x^ + «-^(Vxn 



(47) 




(48) 
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from which it is easy to see that phantom matter eventuaUy comes to dominate the 
universe. Since the Ricci scalar is given by 



there is a future curvature singularity ai t = wtm/i^ + w)- the so-called hig rip. 

It is, however, simple to construct models of phantom energy in which a future 
singularity is avoided. With a suitable choice of potential, we may construct a simple 
scalar field model that exhibits a period of time in which the expansion proceeds with 
w < — 1 and yet settles back to w > — 1 at even later times, thus sidestepping the 
predictions oiw — constant models, and avoiding the big rip. 

However, the background cosmology aside, it is relatively straightforward to see 
that fluctuations in a phantom fleld have a negative energy, and that it may therefore 
be possible for the vacuum to decay into a collection of positive-energy and negative- 
energy particles. Thus, if this instability is present on a timescale shorter than the age 
of the universe, then such a phantom component cannot be a viable candidate for dark 
energy. Such a question has been studied in some detail by several authors [86, 87]. 
For example, in [86], the authors considered a specific toy model in which the null 
dominant energy condition is violated, but in which the cosmology is well- behaved and 
the theory is stable to small perturbations. In many ways, this model has been given 
the maximum opportunity to avoid a problem with instabilities. Nevertheless, what is 
relevant is a field theory calculation of the decay rate of phantom particles into gravitons. 
Clearly, this decay rate would be infinite if the phantom theory was fundamental, valid 
up to arbitrarily high momenta, and would render the theory useless as a dark energy 
candidate. Considering instead the phantom theory to be an effective theory valid 
below a scale A, and including in the Lagrangian operators of all possible dimensions, 
suppressed by suitable powers of the cutoff scale, it was found that such higher order 
operators, even though they may be very high order, can lead to unacceptably short 
lifetimes for phantom particles unless the cutoff scale is less than 100 MeV. This is, of 
course, difficult to reconcile with current accelerator experiments. 

For this, and other reasons, it has proven extremely difficult to construct sensible 
models for which an equation of state w < — 1 is realized within General Relativity 
with minimally-coupled sources, although some progress has been made with ghost 
condensate models [88]. It is, however, important to emphasize that well-behaved 
energy-momentum sources can nevertheless yield an effective equation of state parameter 
that is less than —1 if the dark energy is non- minimally coupled [82], or if the underlying 
gravitational theory is not GR [89, 90, 91]. 

5.2. K-Essence 

In addition to considering the behavior of general potentials for scalar fields, it is also 
possible to think about the effects of modifying the kinetic term, yielding a class of 
models called K-essence [92, 93, 94, 95, 96]. 




(49) 
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To simplify notation in doing this, let us define X = —^g'^'^dij,(f)d^(f). Thus, a general 
action for a scalar field takes the form. 

J d^xyf^L{4>, X) , (50) 

with L an unspecified function, and the canonical scalar field action corresponding to 
the special case L{(j),X) = X — V{(j)). While not completely general, an interesting 
subset of these actions are those with Lagrangian of the form 

L{(f),X) = K{X)-V{c^) , (51) 

with K{X) a positive semi-dcfinitc kinetic function, to be specified. 

If we are interested in background cosmological solutions, with the FRW ansatz 
for the metric, we must assume a homogeneous form for (6, yielding X = ^0^. The 
energy-momentum tensor for this field is straightforward to calculate and yields the 
usual perfect fluid form with pressure p and energy density p given by 

p = C = K{X) - V{<l>) , (52) 

p = [2XKx{X) - K{X)] + V{<f>) , (53) 

where a subscript , X denotes a partial derivative with respect to X Then, deflning 
W(j) = p/p one obtains 

[2XK,x{X) - K{X)] + V{ct>) ■ ' 
Requiring that the energy density of the theory satisfy p > yields 

g{X)-2Xg,x{X)<'^ , (55) 

and a further necessary condition that the theory be stable is that the speed of sound 
of be positive [96] . This yields 

- dp p,x K^X) + 2XK,xx{X) ^ • ^ ^ 

Solving the equations of motion in concert with the resulting Einstein equation, it 
has been shown [92, 93, 94, 95, 96] that tracking behavior can be obtained. This means 
that for a wide range of initial conditions, the energy density of the field naturally 
evolves so as to track the energy density in matter, providing some insight into why 
dark energy domination began only recently in cosmic history. However, as in all rolling 
scalar models, some fine-tuning remains, since one must ensure the right amount of dark 
energy density today. Further, just as with minimal models, it is possible [97] to obtain 
phantom-like behavior, but once again, this comes at the risk of serious instabilities. 

Most interestingly though, the varying speed of sound (compared to = 1 for 
canonical models) opens the window to possible observational signatures of these models 
(which can otherwise be difficult to distinguish from canonical ones [98]) in observations 
of large scale structure in the universe. 
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5.3. Coupled Dark Energy Models 

Our cosmological models seem to require both dark matter and dark energy. A logical 
possibility is that these dark sectors interact with each other or with the normal 
matter [99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109] and, indeed, a number 
of models have been proposed that exploit this possibility to address, for example, the 
coincidence problem. 

There are, of course, many ways in which couplings may arise, and it is not possible 
to provide a comprehensive treatment here. Instead, let us mention an interesting 
subclass of models, defined by the action 



lM|i?-i(V0r-V^(</,) 



+ Sj5j[e2"^W^?,,,vl/j] , (57) 



where g^j^i^ is the Einstein frame metric, is a scalar field which acts as dark energy, and 
are the matter fields. The functions aj{(/)) are couplings to the j*'' matter sector. 
This general action encapsulates many models studied in the literature [110, 111, 
99, 112, 113, 114, 115, 116]. The field equations are 

+ J2e"'^^'^H{Pi+Pi)ui,Ui,+Pig,.] , (58) 
j 

V,W<I> - y (0) = J] ay0)e^«J W(pj - 3pj) , (59) 
j 

where we have treated the matter field(s) in the j*'* sector as a fiuid with density pj 
and pressure pj as measured in the frame e^^'g^,,, and with 4- velocity u-^^j, normalized 
according to g^^u-^^u^^ = — 1. 

An interesting subset of these models is that in which only the dark matter is 
coupled to the dark energy, with the baryonic sector uncoupled. This corresponds to 
setting Q;fe(0) = and q;c(0) = q;(0). Defining pj = e^"jpj gives 

Mp'G^. = V^(/)V,(/) - ^^;,.(V(/.)' - V{(t))g^,, 

+ e^'^'^ p,u,^u,, , (60) 

and V^V^0 — V^^{(i)) = 0, where we have defined an effective potential by Kflf(0) = 
y(0) + e^^^-^^pc ■ The fluid obeys V^(pcO = 0, and <V^< = -{g^''' + «)V^q;. 

There are a number of interesting effects in coupled models, including a surprising 
problem - the adiabatic instability [117, 118, 119, 120, 121]. In the regime in which 
the scalar field adiabatically tracks the minimum of the effective potential, a careful 
study of linear perturbation theory implies that, within a certain regime of parameters, 
exponential growth can occur. 

In general, one can constrain such coupled models using a combination of 
cosmological datasets. In [122] the datasets used included the measurements of the 
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Figure 6. Joint 68% (dark sliaded) and 95% (light shaded) constraints using combined 
WMAP CMB, SDSS matter power spectrum, SDSS and 2dFGRS baryon acoustic 
oscillation, 'union' type la supernovae datasets and HST Ha prior for the power law 
potential (blue) and exponential potential (red) for the the fractional matter density, 
r^M and the effective scalar equation of state, w^, (left panel), and the coupling, C 
(right panel). From [122]. 



CMB temperature and polarization power spectrum from the WMAP 5-year data 
release [123, 39], the 'union' set of supernovae compiled by the Supernovae Cosmology 
Project [37], the matter power spectrum of Luminous Red Galaxies (LRG) as measured 
by the Sloan Digital Sky Survey (SDSS) survey [41, 124] and a Gaussian prior on 
the Hubble constant today, Hq = 72 ± 8, using the Hubble Space Telescope (HST) 
measurements [125]. The resulting constraint plots are shown in figure (6), expressed 
in terms of the dimensionless coupling variable C, defined by 



with /3 = >/2j3. 

6. Modified Gravity 

As we have seen in the previous sections, within GR, a period of accelerated expansion 
can only be due to a dominant component of the energy density with a negative 
equation of state, the most common example being a static, or slowly rolling, scalar 
field (commonly referred to as dark energy or quintessence models). There is, of course, 
a natural alternative to be considered, namely the possibility that the dynamical rules 
determining how the geometry of spacetime responds to regular matter sources may 
differ from those given by GR. In other words, we may consider modifying GR in the 
low curvature regime, such as to admit self-accelerating solutions in the presence of 



negligible matter [126, 127, 128, 129, 130, 110, 131, 132, 133, 134, 135, 136, 137, 138, 
139, 140, 141, 142, 143, 144, 145]. 



C(0) 




(61) 
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While General Relativity has been tested to high precision in laboratory, in the 
solar system and astrophysically, for example with binary pulsars [146], there are far 
fewer tests on cosmological scales. It is therefore natural to explore whether the laws 
of gravity could be modified on the larger scales and whether such modifications could 
be responsible for the late time acceleration. Obviously though, just as with dark 
energy approaches, models of modified gravity need to explain why, or assume a reason 
that the vacuum energy would be zero. With respect to the cosmological constant 
approach however, dark energy and modified gravity models have the advantage of 
adding dynamics to the theory, allowing for the possibility of alleviating the coincidence 
problem and increasing the experimental testability of the models. 

When modifying gravity, one needs to ensure that GR is recovered at early times 
and on short scales, in regimes in which it has been well tested. Moreover, a viable 
model must agree with the precision measurements of the expansion history of the 
universe, of the CMB and of the growth of structure. As we shall see though, on 
theoretical grounds there are obstacles to constructing infrared modifications of GR 
that will escape detection on smaller scales. In a weak-field expansion around fiat 
spacetime, GR describes the propagation of massless spin-2 gravitons coupled to the 
energy-momentum tensor T^,^, including that of the gravitons themselves. But such 
a theory is essentially unique; it has long been known that we can start with a field 
theory describing a transverse-tracclcss symmetric two-index tensor propagating in 
Minkowski space, couple it to the energy-momentum tensor of itself and other fields, 
and demonstrate iteratively that the background metric disappears, leaving instead a 
curved metric obeying Einstein's equation. 

It is therefore generally believed that infrared modifications of GR will necessarily 
involve the introduction of new degrees of freedom. This can be a good thing, in that 
the dynamics of these new degrees of freedom may be precisely what one needs to drive 
the accelerated expansion of the universe. However there are often problems associated 
with the extra dynamics and the extra forces mediated by these degrees of freedom. 
The problems may be of several difi^erent kinds. First, there is the possibility that, 
along with the desired deviations from GR on cosmological scales, one may also find 
deviations on smaller scales; e.g. in the solar system, where GR is well tested. Second 
is the possibility that the newly-activated degrees of freedom may be badly behaved 
in one way or another; either having the wrong sign kinetic terms {ghosts), and hence 
being unstable, or leading to superluminal propagation. More generally, theories that 
contain higher than first order derivatives in the Lagrangian risk being plagued by the 
instability described by the Ostrogradski theorem [147, 148], i.e. the Hamiltonian of the 
system may be unbounded from below. 

The previous arguments put surprisingly restrictive constraints on models of 
modified gravity, and this has motivated some authors to explore models that modify 
GR without introducing new degrees of freedom [129, 149, 150]. 

In this section we will give an overview of some of the main candidate models of 
modified gravity, focusing on the general challenges that such theories face and the 
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peculiar signatures to which they can lead in cosmology. 

There are several ways in which one can imagine of modifying the laws of gravity on 
the largest scales. A first approach could perhaps focus only on the phenomenology of 
cosmology, in which case one postulates specific modifications of the equations describing 
the dynamics of the universe such as to achieve the acceleration. An example is the 
Cardassian model [128], in which the Priedmann equation is modified such that the 
dependence of the Hubble parameter on matter sources is different than the standard 
one. Although models such as this can yield cosmic acceleration, it is desirable to begin 
with a theory described by a generally covariant action, from which we may derive not 
only a modified Friedmann equation, but also the solution describing any other physical 
setting we may wish to know about. 

One broad class of modified gravity models that have been studied in detail contains 
those theories which involve covariant modifications of the 4-dimensional action by 
adding higher order curvature invariants to the Einstein-Hilbert action. The simplest 
examples in this class are the so called f{R) theories. It is instructive to go over some of 
the details of f{R) models since they are representative of generahzed modified theories 
of gravity. We will do so in the following subsection, while more details about these 
models can be found in the recent review by Faraoni et al. [151]. 

Finally, modifications of gravity can originate from extra dimensional scenarios, 
in which all matter fields arc confined to a 4-dimensional brane embedded in a higher 
dimensional bulk. The most prominent example is the Dvali-Gabadadze-Porrati (DGP) 
model [127], in which the 4D brane is embedded in a 5D Minkowksi bulk of infinite 
volume, into which gravity may leak, weakening its strength. We dedicate the last 
subsection to this class of models and to some recent proposals that elaborate on them, 
e.g. Degravitation [152] and Cascading Gravity [153, 154]. 

6.1. f{R) theories 

It is a well known fact that the Einstein-Hilbert action is non renormalizable, therefore 
standard GR cannot be properly quantized. However renormalizability can be cured 
by adding to the action higher order terms in curvature invariants [155, 156, 157, 158]. 
Moreover, higher order terms are generally expected in low energy effective actions 
derived from string theory or Grand Unification Theories (GUT). In 1979 Starobinsky 
showed that one can achieve a de Sitter phase in the early universe by adding to the Ricci 
scalar R in the Einstein action a function f{R) oc R^ [159, 160]. More recently, f{R) 
models have been revisited in the context of cosmic acceleration, starting with [130, 110] 
and followed by many others. The underlying idea is to add functions of the Ricci scalar 
which become important at late times, for small values of the curvature. 
Let us focus on the following action 



where f{R) is a general function of the Ricci scalar, it!. The matter Lagrangian, is 




(62) 
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minimally coupled and, therefore, the matter fields, Xij along geodesies of the metric 
Qi^,^. The field equations obtained from varying the action (62) with respect to g^i, are 

(1 + fn) R,u - \g,u {R + f) + {g^uD - V^V,) fn = , (63) 

where we have defined fn = df /OR. We take the energy-momentum tensor to be that 
of a perfect fiuid (5). The metric Qf^^i, is minimally coupled to matter, hence the stress 
tensor and its conservation law will be the same as in standard GR. 

The f{R) term in the gravitational action leads to extra terms in the Einstein 
equations. The constraint equation of GR becomes a third order equation, and the 
evolution equations now become fourth order differential equations (in contrast to the 
second order equations of standard GR). In particular, when taking the metric to be of 
the flat FRW form (1), the Friedmann equation becomes 

H' + {--fR + Hfn=^P (64) 
6 a 3 

and the acceleration equation is 

--fRH' + { + ^ = -^{p + 2,P). (65) 
a bib 

Cosmic acceleration in equations (64) and (65) can arise through the additional 

terms present for f{R) 7^ 0. One can interpret these extra terms as the contribution 

from an effective fluid with an equation of state 



1 2 
^^^ = -3-3 



H^fn -L-Hfn- Ifn] 

(66) 



It is important to notice that equation (66) can be used as a differential equation for 
the function f{R). Indeed, given an expansion history, with a specified matter content, 
one can determine Wes and then solve eq. (66) for f{R). Since the resulting equation for 
f{R) is of second order, one has a family of f{R) models which reproduce the desired 
expansion history. In this sense, there is a degeneracy among models at the level of 
background cosmology, and, as we will discuss shortly, the growth of perturbations will 
play a key role in breaking this degeneracy. 

There exists a complementary, and sometimes conceptually simpler, way in which to 
approach f{R) theories. It is possible to conformally transform the metric and bring the 
gravitational part of the action to the usual Einstein-Hilbert form of standard GR. The 
price one pays for this simplification is a non-minimal coupling between matter fields 
and the transformed metric [101, 104, 161], as well as the appearance of a new scalar 
degree of freedom playing the role of DE and evolving under a potential determined by 
the original form of the function f{R) in (62). When the transformation between the 
two frames is well-defined, the classical results obtained using either description should 
be the same. The frame in which matter particles fall along geodesies of the metric, but 
in which the gravitational part of the action is modified, is referred to as the Jordan 
frame. The frame in which the gravitational part of the action is the same as in standard 
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GR, while the matter may be non-minimally coupled to gravity, is called the Einstein 
frame. 

The emergence of an additional scalar degree of freedom can also be seen directly 
in the Jordan frame. The role of the scalar field is played by f^, dubbed the scalaron 
in [160], and obeys an equation of motion given by the trace of equation (63), which 
can be written as 

□/« = 3 (i? + 2/ - Rfn) - _(p - 3P) ^ ^ , (67) 

which is a second order equation for /^j, with a canonical kinetic term and an effective 
potential Kff(/R)- Any f{R) theory designed to achieve cosmic acceleration must satisfy 
1/ -C -R| and -C 1 at high curvatures to be consistent with our knowledge of the 
high redshift universe. In this limit, the extremum of the effective potential lies at the 
GR value R = K^{p — 3P). Whether this extremum is a minimum or a maximum is 
determined by the second derivative of Ves at the extremum, which is also the squared 
mass of the scalaron 



m 



2 _ d^Ves 1 



1 + /^ p' 



RR 



(68) 



At high curvatures, when \RfRR\ <^ 1 and /r — > 0, 



.2 ^A + fR^. 1 



~ -Tr~ ~ w~ ■ ^^^^ 

It then follows that in order for the scalaron not to be tachyonic, (i.e. to have a positive 
squared mass), one must require fRR > 0. Classically, fRR > is required to keep the 
evolution in the high curvature regime stable against small perturbations [162, 163]. 

The scalaron mediates an attractive "fifth force" , which has a range determined by 
the Compton wavelength 

K^^. (70) 

While Ac is large at current cosmological densities, terrestrial, solar and galactic tests 
are not necessarily violated because the scalaron acquires a larger mass in regions of 
high matter density. This is essentially the Chameleon mechanism of [112, 113]. 

As we mentioned above, alternatively one may map the theory to the Einstein frame; 
following the approaches of [111, 161, 164], one can recast the gravitational action (62) 
into a dynamically equivalent form by applying the conformal transformation 

~9,. = e'-^^'-^g,. , (71) 

such that the function uj{x°') satisfies 

e-'^{l + fR)^l . (72) 

With this choice of u; the action (62) transforms into an action with the usual Einstein- 
Hilbert form for gravity. If we now define the scalar field = 2uj/I3k, where (3 = •\/2/3, 
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the resulting action becomes 
1 
2^ 



-5= At / d''xJ-gR+ I d^xJ-'g 



-lg>^^{V,(t>)^A-V{(t>) 



+ j d'x^ge-^^^'^U[xue-^^'>'~g,J, , (73) 



where the potential V{(j)) is determined entirely by the original form (62) of the action 
and is given by 

and where a tilde is used to indicate Einstein's frame quantities. 

The field equations obtained by varying the action with respect to the metric g/^^, 

are 

G^, = StiGT^, + ^^cpV^cp + ^ir^Va<l>S/^<l>)g^. - V{(P)~g^, (75) 



{f,^ + T^) , (76) 

where T^^ is the energy-momentum tensor of matter fields in the Einstein frame. In the 
remaining terms on the right hand side of (75) we recognize the energy-momentum tensor 
T^^ of the scalar field 0. Although quantities in the Einstein frame are more familiar 
than those in the Jordan frame, there are some crucial distinctions. Most notably, in 
this frame test matter particles do not freely fall along geodesies of the metric ^^j,, since 
the scalar field is also coupled to matter, as can be seen from the conservation equations 
for the scalar field and for the perfect fluid matter, which read 

+ + l^/3(p_3P) , (77) 

'p + 3H{p + P) = -\k134>{p - 3P) , (78) 
where = dV/d(f). 

From this mapping it is clear that f{R) theories are a specific example of the 
broader class of scalar-tensor theories that we discussed in the final subsection of the 
previous section. 

Returning to the Jordan frame formulation, there is a series of conditions that must 
be imposed on the function f{R) in order to avoid instabilities and ensure agreement 
with local tests of gravity. These conditions can be summarized as follows 

(i) Irr > Q ioi R ^ fuR. Classically, this follows from requiring the existence of a 
stable high-curvature regime, such as the matter dominated universe [162]. 

(ii) 1 -I- /k > for all finite R. The most direct interpretation of this condition is that 
the effective Newton constant, Gefr = G/{1 + /r), is not allowed to change sign. 
This prevents the graviton from becoming ghost-hke [165]. 

(iii) fa < 0. Given tight constraints from Big Bang nucleosynthesis (BBN) and the 
CMB, we must require GR to be recovered at early times, i.e. f{R)/R — > and 
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/ij — > as i? — s> oo. Together with > 0, this imphes that fa must be a negative, 
monotonicaUy increasing function of R that asymptotes to zero from below. 

(iv) fjt must be small at recent epochs. This is not required if the only aim is to build a 
model of cosmic acceleration, without trying to satisfy the solar and galactic scale 
constraints. This condition ensures a small difference between the value of the 
scalaron field in the high density galactic center and the value in the outskirts of 
the galaxy [166]. As shown in [113], the difference between these values is effectively 
the potential difference that sources the attractive "fifth force" acting on objects 
in the vicinity of the galaxy. The analysis of [166] suggests that the value of 
today should not exceed 10~^. However, as stressed in [166], this bound assumes 
that in f{R) models galaxy formation proceeds similarly to that in GR. Hence, 
while it is certain that must be small, any specific bound on its value today 
will be unreliable until galaxy formation in f{R) is studied in N-body numerical 
simulations - work that has only recently begun [167, 168]. 

Finally, it was recently shown [169, 170] that even models of f{R) that satisfy all 
of the above conditions might still face problems in the strong gravity regime. These 
problems are due to a curvature singularity R ^ appearing at the non-linear level, at 
a finite field value and energy, such that the scalar degree of freedom can easily access it 
in the presence of matter. As a consequence, f{R) theories appear to be incompatible 
with the existence of static dense compact objects like neutron stars. The problem 
seems to require an ultraviolet completion of f{R) theories, and indeed it was recently 
shown in [171, 172] that the addition of higher order curvature corrections might solve 
the curvature singularity problem, although an unpleasant degree of fine-tuning would 
be required. 

6.2. Generalized theories of gravity 

The f{R) theories studied in the previous subsection are the simplest example of 
generalized theories of gravity. Still, in their simplicity, they retain some of the features 
of more general higher order theories of gravity. Studying the cosmology of these models, 
we have learned how the introduction of higher order terms in curvature invariants can 
lead to instabilities and/or confiicts with local tests of gravity because of extra degrees of 
freedom that come into play. These problems are characteristic of generalized theories of 
gravity, and can be more severe in models that involve functions of curvature invariants 
other than the Ricci scalar. 

There are, of course, any number of terms that one could consider, but for simplicity, 
we might consider those invariants of lowest mass dimension that are also parity- 
conserving 




(79) 
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and actions of the form [142] 

S = J d*x^[R + f{R,P,Q)] + J d^x^CM, (80) 

where P, Q) is a general function describing the deviations from general relativity. 
It is convenient to define 

^^=a^' ^^=ap' ^^ = aQ' ^^^^ 

in terms of which the equations of motion are 

+ Jr R^v + 2/p -R"^ Rav + 2/q RaP^n R°'^''v 

+ VaV;3(/pi?"'^) - 4V«V,fl[/Q = SttGT^, . (82) 

It is straightforward to show that actions of the form (80) generically admit a 
maximally-symmetric solution: R — a non-zero constant. However, an equally generic 
feature of such models is that this de Sitter solution is unstable (corresponding to one 
of the new degrees of freedom attaining a local potential maximum) , and in many cases 
there exists another accelerating solution, analogous to those that we found for f{R) 
models. 

What about constraints on these models? It has been shown [173] that solar system 
constraints, of the type wc have described for f{R) models, can be evaded by large classes 
of these more general models. Roughly speaking, this is because the Schwarzschild 
solution, which governs the solar system, has vanishing R and P, but non- vanishing Q. 

Far more serious though is the issue of ghosts and superluminal propagation. It 
has been shown [174, 175] that a necessary but not sufficient condition that the action 
be ghost-free is that b — —4:C, so that there are no fourth derivatives in the linearized 
field equations. What remained was the possibility that the second derivatives might 
have the wrong signs, and also might allow superluminal propagation at some time 
in a particular cosmological background. It has recently been shown that around a 
FRW background with matter, the theories are ghost-free, but contain supcrluminally 
propagating scalar or tensor modes over a wide range of parameter space [176, 177]. 
However, more generally it has been argued [148], that only modifications of the action 
that are functions of the Ricci scalar and the Gauss-Bonnet term can avoid the instability 
of the Ostrogradski's theorem. Thus, such general actions do not seem to be promising 
candidates for modified gravity theories explaining cosmic acceleration. 



6.3. Modifying gravity without extra degrees of freedom 

The challenges that higher order theories of gravity theories face are related to the 
dynamics of, and forces mediated by, the extra degree (s) of freedom introduced by 
the modifications. It is therefore interesting to investigate whether it is possible to 
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achieve cosmic acceleration with models for which the constraint equations of gravity 
arc modified while the graviton is still the only propagating degree of freedom [129, 149]. 
Since such models are a little outside the scope of this review, we will not comment on 
them further, other than to point out that linear perturbation theory seems to show a 
possible conflict with observations of the formation of structure [129]. 



6.4- The DGP model and its derivatives 

So far we have considered models of modified gravity based on a covariant four- 
dimensional action. In this subsection we shall review some of the proposals in the 
context of extra dimensions. We focus on brane-world models that modify gravity at 
low energies; i.e. models that can be used to achieve cosmic acceleration. 

The most prominent proposal in this context is the well known Dvali-Gabadadze- 
Porrati (DGP) [127] model, the cosmology of which was described in [126]. In this 
model, the matter fields are constrained to lie on a brane embedded in a fiat 5D 
Minkowski bulk of infinite volume, into which the graviton can leak. In the DGP model 
gravity appears four dimensional at short distances, while it is modified on distances 
large compared to a crossover scale Tc through the evaporation of the graviton into the 
unseen fifth dimension. 

The 5D action describing the model is 



—M^ f d^xyf^R ^—Ml I d^x y/-g(^^ 

.67r J 16n J 



^ Ml 



167r 

where M and R represent, respectively, the five dimensional metric, Planck mass and 
Ricci scalar, while g^'^\ Mp and it!^^) are their four dimensional counterparts. Finally, 
is the AD matter Lagrangian density. The crossover scale Tc is given by the ratio of 
the observed 4D Planck mass to the 5D one 

^ - (84) 



(83) 



2M3 ■ 

This represents the distance over which metric fiuctuations propagating on the brane 
dissipate into the bulk, weakening gravity. As a consequence, the 4D effective 
gravitational potential will behave as 

$oc|^"' . (85) 

Interestingly, significant departures from standard gravity persist down to scales 
smaller than r^; for a mass source of Schwarzschild radius r^,, modifications will be 
relevant down to scales equal to the Vainshtein radius r* = {rgrlY^^ [178, 179]. 

In order to address the phenomenon of cosmic acceleration it is necessary to choose 
Tc ~ Hq^. The Priedmann equation describing the expansion of a 4D brane with spatial 
curvature k reads 



k 1 k 87rG' 



i/^ + - - + p. (86) 



rc V 
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where e = ±1 correspond to the two different branches of the modeL The e = +1 choice 
corresponds to the self-accelerated branch, sDGP; this model naturally leads to the late 
time acceleration and has been recently tested against SNcIa, BAO and CMB data [180]. 
It was found that the flat sDGP {k = 0) model, is a poor fit to the data and that the 
ACDM model gives a significantly better fit [89, 181, 182, 183, 184]. Moreover, sDGP 
is plagued by the presence of a ghost in the gravitational background sector, which is a 
serious problem. 

The choice e — —1 corresponds to the normal branch, nDGP, which is obtained 
via a different choice of embedding of the brane into the bulk. It requires the addition 

of a cosmological constant to achieve cosmic acceleration, but it has the advantage of 
being ghost free. In this model, the cosmological constant is screened by 5D effects, and 
therefore it is the simplest example of degravitation, which we will mention shortly. 
However the screening effect is too small to solve the fine tuning problem of the 
cosmological constant. Another interesting feature of nDGP is that it displays an 
effective phantom behavior [89, 181, 182, 183]. A comparison of the nDGP model 
with a combination of cosmological datasets [89, 181, 182, 183, 184], found that closed 
models were a good fit to the data. 

In [152] it was pointed out that massive-graviton theories, such as DGP, could lead 
to degravitation of the vacuum energy, offering a dynamical solution to the cosmological 
constant problem. Degravitation indeed allows for a large cosmological constant, but 
suppresses its backreaction by making gravity sufficiently weak on large scales. The 
degravitation of the vacuum energy is an intriguing solution to the fine tuning problem 
of the cosmological constant, and can be achieved if Newton's constant is promoted to 
a derivative operator [152] 

G-^ {L^n) G^, = 87tT^,. (87) 

The gravity described by the above Einstein equations will behave as a high-pass filter 
with characteristic scale L, such that G — ^ G^^^ for □) oo, and G ^ 
for □) 0. Equation (87) should be considered only as a phenomenological 

description of the degravtitaion model, the details of which can be found in [152]. 

Although DGP is effectively a massive graviton theory, the degraviation 
phenomenon cannot be achieved in the context of standard DGP because the weakening 
of gravity is not sufficiently effective. A very promising way of realizing degravitation 
is through the so called Cascading Gravity model [153, 154, 185, 186]. In this model, 
DGP is extended to higher dimensions, in which our visible 3-brane is embedded within 
a succession of higher-dimensional branes, each with their own induced gravity terms, 
embedded in a flat D-dimensional bulk. The model is free of ghost-like instabilities 
(common to previous attempts at a higher-codimenion DGP) and the 4D propagator is 
well behaved when calculated on the brane. Moreover, this model exhibits degravitation 
with a steep enough weakening of gravity - the gravitational force law cascades from 
a 1/r^ (4D) behavior, to 1/r^, (5D) to 1/r^, (6D) behavior, etc. as one probes larger 
distances on the 3-brane. 
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These models are works in progress, but are particularly interesting examples of 
modified gravity models in which both cosmic acceleration and the cosmological constant 
problem may be addressed. 

7. Distinguishing between models 

The ACDM model is currently the best fit to the available cosmological data, and 
passes all local tests of gravity [146]. As we have already discussed, it is however 
important to explore the whole space of explanations that fit the data comparably well 
and have testable features. In the previous sections we saw that, in general, imposing 
cosmological and local viability of models of modified gravity, as well as generalized DE, 
requires a degree of fine-tuning. After this tuning, the predictions of these models 
for the expansion history of the universe are often identical to that of the ACDM 
model [152, 166, 187, 188, 189, 190]. 

However, this degeneracy is typically broken at the level of cosmological structure 
formation; indeed, models of modified gravity that closely mimic the cosmological 
constant at the background level can nevertheless yield significantly different predictions 
for the growth of structure. In recent times it has become increasingly evident, that 
the large scale structure of the universe offers therefore a promising testing ground 
for GR. To this extent, it is important to explore to what degree one can use LSS to 
distinguish among all the candidate models of cosmic acceleration and, indeed, whether 
model-independent approaches are possible [191, 192, 193, 194]. 

Scalar metric perturbations in the Newtonian gauge are described by two potentials, 
^(x, t) and t), defined via the perturbed FRW line element 

ds'^ = - (1 + 2^(f , t)) df + a^it) (1 - 2$(,f, t)) dx^ . (88) 

Thus, these potentials correspond, respectively, to perturbations in the time-time and 
space-space components of the metric tensor, representing the strength of gravity and 
the spatial curvature. In the ACDM model, the potentials are equal during the epoch 
of structure formation, and their time dependence is set by the same scale-independent 
linear growth function that describes the growth of matter density perturbations. This, 
it turns out, is a very pecuhar feature of the ACDM model, and of models of dark 
energy with negligible shear and clustering. However, it no longer holds in theories of 
modified gravity or models of coupled dark energy, where one can have scale-dependent 
growth patterns. In these theories, the two Newtonian potentials generically differ, 
and their dependence on matter perturbations can be different. As a consequence, 
modifications of gravity generally introduce a time- and scale-dependent slip between 
the Newtonian potentials, as well as a time- and scale-dependent effective Newton's 
constant describing the clustering of dark matter. These modifications are expected to 
leave distinct imprints on the large scale structure of the universe, which may help to 
break the degeneracy that characterizes models of cosmic acceleration at the background 
level. This explains the recent growing interest in studying the potential of current and 
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upcoming surveys to detect and constrain any departure from GR in the growth of 
structure [195, 196, 197, 198, 199, 200, 201, 202]. 

7.1. Modified Growth of Structure 

It is useful to review the evolution of linear scalar perturbations in general models of 
modified gravity or coupled dark energy, as it will allow us to describe in detail some of 
the recent efforts to constrain departures from GR using the growth of structure. 

We focus on the evolution of linear matter and metric perturbations in a general 
metric theory of gravity. Assuming that the background evolution is correctly described 
by the flat FRW metric, we consider scalar perturbations in Newtonian gauge, as deflned 
in (88). Representing all perturbed quantities in Fourier space, we use the standard 
notation for the energy momentum tensor of the matter fields, which to first order in 
the perturbations, assumes the form 

T^^-{p + P)vi, (89) 
Tj = {P + SP)6] + 7ii, (90) 

where 6 = Sp/p is the density contrast, v the velocity field, 6P the pressure perturbation 
and TTj denotes the traceless component of the energy-momentum tensor. Finally, we 
define the anisotropic stress a via (p + P)a = —{k''kj — If^pTr]. 

In GR, the linearized Einstein equations provide two independent equations relating 
the metric potentials and matter perturbations, the Poisson and anisotropy equations, 
respectively: 

.2 



*^* = -5^pA. (91) 

*'(4>-*) = ^(p + ^'y, (92) 

where pA = p8 + 3^(p + P)v is the comoving density perturbation. In the ACDM 
and minimally coupled quintessence models, the anisotropic stress is negligible at times 
relevant for structure formation, and we have ^ = $. 

In models of modified gravity, as well as in more exotic models of dark energy, the 
relation between the two Newtonian potentials, and between the potentials and matter 
perturbations, can be different [203, 204, 205] . The effect of these modifications can be 
encoded in two functions of time and space, a rescaling of the Newton's constant p(a, k) 
and a gravitational slip 7(0, k), defined via 

/^'* = -2^M«,%A (93) 

| = 7(a,A;), (94) 

where /i(a, k) and 7(0, k) are time- and scale-dependent functions. Note that p is defined 
via the Poisson equation (93) written in terms of the perturbation to the time-time 
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component of the metric. This choice is natural, as it is ^ that enters the evolution 
equation for CDM density perturbations on sub-horizon scales: 

5 + 2H5+ = . (95) 

With the set of equations (93), (94) and (95), completed by an expression for fi and 7, 
one can evolve the linear perturbations forward in time and extract predictions for the 
growth of structure. 

7.1.1. An example It is instructive to look at the growth of structure in a specific model 
of modified gravity. For this purpose, we will briefly review some results for the evolution 
of linear perturbations in f{R) theories. The growth of structure in these models has 
been studied in [188, 205, 206], where it was found that the perturbations grow with a 
characteristic scale-dependent pattern. The Compton wavelength of the scalaron (70) 
introduces a scale which separates two regimes of sub-horizon gravitational dynamics 
during which gravity behaves differently. On scales A ^ Ac, the scalaron is massive 
and the "flfth force" it mediates is exponentially suppressed; thus deviations from GR 
are negligible. However, on scales inside the Compton radius, the scalaron is light and 
deviations are significant. The relations between $ and and the relation between 
them and the matter density contrast, will be different below the Compton scale and 
that affects the growth rate of structure. Indeed, on sub-horizon scales {k S> aH), as 
an effect of the modifications to the Einstein-Hilbert action, the Poisson and anisotropy 
equations read 



1 k'^ 

2 a'^ 



(96) 



_ 3 ; ^\ 3 

2- 

il + fR)i<t>-^) = fRR6R, (97) 

where we have neglected any shear from matter fields. From the sub-horizon limit 
of equations (96), (97), we can easily read off the resulting rescaling of the Newton's 
constant and the gravitational slip as 

a^{l + fR) + SfRRk^ y + XI 
,^ a\l + fR)+2fRRk^ 3A^ + 2X1 .qq^ 
- a^{l + fn) + AfRRk^ - 3P + 4A^ ' ^^^^ 

(where we have neglected a prefactor (1 -|- /h)"^ multiplying the expression for since 
it is very close to unity in models that satisfy local tests of grvaity [166]). In particular, 
on scales below Ac, the modifications introduce a slip between the Newtonian potentials, 
leading to the low-curvature regime ^ ^ 2$. The rate of growth of structure depends 
on the balance between the fifth force and the background acceleration. For modes that 
cross below Ac during matter domination, the effect of the modifications is maximized 
as the potentials grow in the absence of background acceleration. When, however, the 
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Figure 7. The function probed by the cross-correlation of large scale structure 
with ISW, 6 ■ + '^)/dz, as a function of scale k and redshift z. The left panel 
corresponds to an f{R) model with a ACDM expansion history and f'^ = —10^^ 
and shows a characteristic scale-dependent pattern. The right panel corresponds to 
ACDM. The dashed line crossing through the left panel corresponds to A = Ac (70), 
i.e. it corresponds to the Compton wavelength of the scalaron Ac (70). For scales 
A < Ac, during matter domination, one can clearly notice the effect of the "fifth force" 
which suppresses the cross correlation and can actually make the correlation negative. 
Therefore, a negative cross correlation signal at early redshifts (corresponding to matter 
era) , is a signature of / (i?) . The acceleration of the background will eventually contrast 
the "fifth force" and lead to a positive cross-correlation. 

background starts accelerating, the potentials begin to decay but at a lesser rate than 
in the ACDM model. Therefore, a characteristic signature of f{R) theories would be a 
non-zero ISW effect during the matter era [205, 188]. The ISW constitutes just a small 
fraction of the overall CMB power on the larger scales, where cosmic variance dominates, 
and therefore it is very hard to extract. The way to proceed is to cross-correlate the 
CMB with galaxies, since one expects a correlation between the distribution of galaxies 
(i.e. DM potential wells) and the ISW signal [48, 207, 208]. In f{R) theories such a 
correlation would be negative, since the potentials would be growing. In Fig. (7) we 
show such cross correlation for an f{R) model, while some constraints using existing 
ISW data can be found in [209]. One should note, however, that the ISW measurements 
are plagued by large statistical errors (because the ISW effect is only a part of the total 
CMB anisotropy) and cannot provide the percent level accuracy needed to test viable 
f{R) models, (i.e. models that satisfy local tests of gravity and for which the ISW 
contribution is suppressed to the percent level). Weak lensing studies, on the other 
hand, are counted on to eventually provide highly accurate 3D maps of the sum of the 
gravitational potentials, $ -f 

The features found in f{R) models are common to the broader class of scalar-tensor 
theories, to which f{R) belongs, as well as to models of coupled quintessence [210, 105]. 
These models can be described by a Chameleon-type action (57), representative of 
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coupled dark energy models in the Jordan frame and of scalar-tensor theories in the 
Einstein frame. Focusing, as before, on dark matter, indicating its coupling simply 
with a (i.e. dropping the index i), and studying scalar linear perturbations around a 
flat FRW metric, on sub-horizon scales we can easily find [188, 205] that the effective 
Newton's constant and the gravitational slip are 



1 + (1 + 2Mj,a%) ^ 



2 



/.(a, k) = I. ° " (100) 

1 + (1 - 2M2ay J 



7(«, k) = — ) , (101) 

l + (l + 2M>f,)^ 

which arc exactly of the form (98)-(99). (In (100) we have neglected e prefactor 
g-2a(<^) ^ 2). Here m is the mass of the scalar field and a ,^ = da/dcf). 

To summarize, in models of modified gravity and coupled dark energy the growth 
pattern is in general scale-dependent and the time evolution of the Newtonian potentials 
is different than that in ACDM. The introduction of an effective Newton's constant and 
a slip between the potentials is a general feature of the modifications, independent of 
the background expansion. Therefore, even for models that at the background level 
closely mimic the cosmological constant, the predictions for the growth of structure can 
be significantly different. The degeneracy characterizing models of cosmic acceleration 
at the background level can henceforth be broken at the level of large scale structure. 
Characteristic imprints of such modifications will be a different clustering rate for DM, a 
modified ISW signal as well as a different Weak Lensing shear (WL). Moreover, since the 
potentials $ and ^ will in general differ, one expects that the function measured by ISW 
and WL, i.e. ($-1-^'), will be different than the one to which peculiar velocities respond, 
i.e. ^. All these features allow potentially powerful ways to distinguish between models 
of cosmic acceleration, which we review in the next subsection. 

Before concluding, let us notice that, as we saw, models of generalized dark energy 
can induce a pattern of growth analogous to the one found in models of modified gravity. 
Therefore, the degeneracy between dynamical dark energy and modified gravity models 
appears to persist at the level of LSS. Indeed, it seems always possible to introduce a 
generic dark fiuid with properties such as to mimic the modifications of gravity. This has 
been discussed recently by several groups [211, 212, 213] and the authors of [213] have 
noted that the breaking of such degeneracy would require the independent measurements 
of three observables, (e.g. ^, $ -|- ^ and 5). 

7.2. Constraining departures from standard GR 

From the preceding discussion it is clear that LSS offers a potentially powerful way of 
constraining models of cosmic acceleration. The task, however, is not trivial, and there 
are different ways in which one can proceed. A first approach could perhaps consist of 
assuming a functional form for the slip between the potentials and the effective Newton's 
constant, motivated by a certain class of theories. Some representative parametrizations 
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have been introduced in [197, 212, 214, 215]. From the previous examples we can easily 
infer that a quite generic parametrization would be 

where the parameters Af have dimensions of length squared, while the f3i represent 
dimensionlcss couplings. Finally, from a scalar-tensor point of view, the parameter 
s encodes the time-dependence of the scalar field mass. The expressions in (102) 
coincide with the scale-dependent parametrization introduced in [197]. In [201], the 
authors performed a Fisher matrix analysis to forecast the constraints on the parameters 
{Aj,/3j,s} of (102) from a combination of power spectra from galaxies. Weak Lensing 
and the Integrated Sachs- Wolfe effect on the CMB. Such an analysis reveals the extent 
to which one can constrain these theories and also allows one to reconstruct the shape of 
the slip and effective Newton's constant based on the chosen form. In [201], it was found 
that these data, even at the linear level, are quite powerful in constraining the modified 
growth parameters. The results, of course, depend on the choice of the parametrization. 
However, they are good indicators of the power of current and upcoming surveys to 
constrain departures from GR. 

An alternative approach, which can work for certain estimators of the slip, is a 
direct reconstruction from data. In [199, 200], it was proposed to consider the ratio 
of the peculiar velocity-galaxy correlation with the Weak Lensing-galaxy correlation. 
Specifically, the authors propose the following estimator for the gravitational slip 

(104) 

where Ps/^(^^l-<s>)g and P\r2^g are the cross-power spectrum between the galaxy number 
overdensity and, respectively, weak lensing and the potential ^ (which is related to 
peculiar velocities via the continuity equation). 

In such a ratio, the dependence on the galaxy bias cancels out. Such a ratio, 
if appropriately constructed, would directly probe any difference between $ and ^. 
This is a more direct and model-independent way of testing GR with the growth of 
structure; however its power will depend on how well future experiments will be able 
to measure peculiar velocities. From the analysis in [199, 200], it appears however that 
the ratio measured by future surveys might have strong discriminating power for some 
dark universe scenarios. 

Finally, another, non-parametric approach, consists of performing a Principal 
Component Analysis (PCA) to determine the eigenmodes of the slip and the Newton 
constant, that can be constrained by data, in the same way as has been done for the dark 
energy equation of state [216, 217, 218]. This method allows one to compare different 
experiments and their combinations, according to the relative gain in information about 
the functions. PCA can also point to the "sweet spots" in redshift and scale where data 
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is most sensitive to variations in the slip and Newton's constant, which can be a useful 
guide for designing future observing strategies. The PCA method does not allow one to 
reconstruct the shape of the functions from data. However, one can still reproduce the 
errors on parameters of any parametrization from the eigenvectors and eigenvalues found 
using PCA [217]. Thus, it is a promising method, but computationally challenging. 

7.2.1. The observables As we mentioned above, in the ACDM model the sub-horizon 
evolution of gravitational potentials and the matter density fluctuations are described 
by a single function of time - the scale- independent growth factor g{a). In models 
of modified gravity, on the other hand, the dynamics of perturbations can be richer 
and, generically, the evolution of $, and S (the matter density contrast), will be 
described by different functions of scale and time. Therefore, we expect that different 
observables will be described by different functions, and by combining different types 
of measurements, one can try to reconstruct these functions, or at least put a limit on 
how different they can be. In what follows, we shall give a brief review of the relation 
between the different types of observables and the gravitational potentials they probe. 
For a more thorough overview of the various ways of looking for modifications in the 
growth of perturbations we refer the reader to [213]. 

Galaxy Counts (GC) probe the distribution and growth of matter inhomogeneities. 
However, to extract the matter power spectrum, one needs to account for the bias, which 
typically depends on the type of galaxies and can be both time- and scale-dependent. On 
large scales, where non-linear effects are unimportant, one can use a scale- independent 
bias factor to relate galaxy counts to the total matter distribution. This relation 
becomes increasingly complicated and scale-dependent as one considers smaller and 
smaller scales. In principle, the bias parameters can be determined from higher order 
correlation functions [219, 220, 221]. On sub-horizon linear scales, the evolution of 
matter density contrast is determined by equation (95). Hence, measurements of GC 
over multiple redshifts can provide an estimate of ^ as a function of space and time, 
up to a bias factor. A more direct probe of the potential ^, would be a measurement 
of peculiar velocities, which follow the gradients of ^. Such measurements would be 
independent of uncertainties associated with modeling the bias. 

Peculiar velocity surveys typically use redshift-independent distance indicators to 
separate the Hubble flow from the local flow, and nearby SNela are therefore good 
candidates; a number of surveys, like the 6dFGS [222] and the 2MRS [223], use 
galaxies. An interesting alternative is offered by the kinetic Sunyaev-Zel'dovich effect 
in clusters [224], that arises from the inverse Compton scattering of CMB photons off 
high-energy electrons in the clusters. This effect provides a useful way of measuring 
the bulk motion of electrons in clusters, hence the peculiar velocity of clusters, but it 
is limited by low signal-to-noise ratio. Current measurements of peculiar velocities are 
limited in accuracy, and at this point it is not clear how to forecast the accuracy of 
future observations. Therefore we did not include them in our observables, even though 
they are a potentially powerful probe [225] . 
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In contrast to galaxy counts and peculiar velocities, which respond to one of the 
metric potentials, namely ^, Weak Lensing of distant light sources by intervening 
structure is determined by spatial gradients of the sum ($ + ^). Hence, measurements of 
the weak lensing shear distribution over multiple redshift bins can provide an estimate 
of the space and time variation of the sum of the two potentials. In the ACDM 
and minimally coupled models of dark energy, the two metric potentials coincide, and 
therefore WL probes essentially the same growth function that controls the evolution 
of galaxy clustering and peculiar velocities. In models of modified gravity, however, 
there could be a difference between the potentials, corresponding to an effective shear 
component. Measurements of the Integrated Sachs-Wolfe effect in the CMB probe the 
time dependence of the sum of the potentials: d ($ + ^) /dt. 

By combining multiple redshift information on GC, WL and CMB, and their cross- 
correlations, one can constrain the differences between the metric potentials and the 
space-time variation of the effective Newton constant defined in the previous section. 
Ideally, experimentahsts would measure all possible cross-correlations, between all 
possible pairs of observables, in order to maximize the amount of information available to 
us. In practice, however, it can be difficult to obtain these cross-correlations, since their 
measurements require that each of the individual fields (CMB, GC, WL) be measured on 
the same patch of sky. This will be addressed with near and distant future tomographic 
large scale structure surveys (such as DES [31], LSST [32] and PAN-STARRS [226]). 
Even with conservative assumptions about the data, (i.e. considering only linear scales), 
it is hoped [201], that DES will produce non-trivial constraints on modified growth, 
and that LSST will do even better (with pa 10% relative errors in the parameters). 

Prom the discussions in the previous section, it is clear that the predictions for 
the growth of structure in models of modified gravity and coupled dark energy can 
significantly differ from those in ACDM or uncoupled quintessence models. Therefore, 
LSS and CMB data offer a valuable testing ground for gravity and cosmological tests are 
expected to play an important role in determining the physics of cosmic acceleration. 
In principle, a combination of galaxy number counts and weak lensing measurements, 
along with CMB and SNe data, will allow a distinction between modified gravity /exotic 
dark energy models and the standard ACDM model of cosmology. However, as we 
have already mentioned, it will be harder to break the degeneracy between models of 
modified gravity and generalized models of dark energy, where the dark fiuid is allowed 
to cluster and carry anisotropic stress. As discussed in [213], the latter task will require 
the independent measurements of at least three observables. 

8. Conclusions 

That the expansion rate of the universe is accelerating is now, just over a decade after 
the first evidence from observations of type la supernovae, a firmly established aspect of 
cosmology. The rapid progress in establishing this fact is a testament to the breathtaking 
convergence of techniques and technology that has emerged in observational cosmology. 
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In turn, cosmic acceleration has introduced new wrinkles into almost every part of 
theoretical cosmology, ranging from the details of structure formation, through the 
CMB and gravitational lensing. 

Beyond the fascinating new problems in theoretical cosmology, cosmic acceleration 
has presented an enormous and, as yet, unmet challenge to fundamental physics. A 
perfectly good fit to all known cosmological data is that cosmic acceleration is driven by 
a cosmological constant. However, accounting for the necessary magnitude of such an 
object seems at least as difficult as attempts to understand why it should be precisely 
zero. Anthropic arguments, fueled by results from the landscape of string theory and the 
idea of eternal inflation, may yield a way to understand such a small number. However, 
it is, at present, far too early to know if this is a sensible outcome of string theory, 
and there are no developed ideas of how such proposals might be tested. This is not 
to say, of course, that they are incorrect; merely that there are immense obstacles to 
establishing such ideas as a tested answer. 

If cosmic acceleration is not due to a cosmological constant, then an exotic 
component of the energy budget of the universe may be the culprit. It is important 
to recognize up front that any such understanding of cosmic acceleration assumes that 
the cosmological constant itself is either precisely zero, or at least negligibly small - 
an unsolved problem in itself. If we make such an assumption, then proposals for a 
dynamical component of the cosmic energy budget, driving the accelerating expansion, 
come under the heading of dark energy. In their simplest forms, ideas of dark energy 
can be thought of as attempts to construct a model for late-time cosmological inflation, 
with the key differences being the extremely low energy scales involved, and the fact 
that this period of acceleration need not end in the way that inflation must. As with 
early universe inflation, a major challenge is flnding a home for the abstract idea within 
a well-motivated and technically natural particle physics model. In fact, in a number 
of ways, this problem is more acute for dark energy models, since they must operate 
in a regime in which we already have a supremely well-tested low-energy effective fleld 
theory of much of the matter content of the universe. As we have described, some 
tentative ideas exist to tackle this problem, but as yet no convincing particle physics 
implementation of dark energy exists, in contrast to the somewhat natural occurrence 
of dark matter candidates in theories of physics beyond the standard model. 

But dark energy is not the only possible solution to the problem of cosmic 
acceleration. Instead we might imagine leaving the material contents of our theories 
untouched, comprised only of the standard model and dark matter, and instead 
revisiting the dynamical relationship between these components and the evolution 
of the background geometry of the universe. In other words, we might consider 
that cosmic acceleration is our first evidence for a modification of General Relativity 
in the far infrared. One might imagine that this approach would hold out the 
possibility of a solution to the cosmological constant problem itself, as well as 
providing an understanding of cosmic acceleration. However, to date this possibility 
remains unrealized, and just as for dark energy, modified gravity approaches to cosmic 
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acceleration must also merely assume that the cosmological constant is either zero, or 
negligibly small. 

Making this assumption, a number of authors have considered how one might 
modify gravity. However, to date all attempts have run into problems, either with 
matching constraints on gravity at small (solar system) scales, or with the existence 
of ghost degrees of freedom, raising the possibihty of problems similar to those we 
mentioned when discussing phantom fields. Nevertheless, if one views these models 
either as a long-wavelength approximation to a better-behaved theory of gravity, or 
as a very low energy effective theory, in need of an ultraviolet completion to cure 
its problems, one may extract a variety of predictions for signatures in upcoming 
cosmological observations. 

Perhaps the biggest question to be answered is binary in nature - is cosmic 
acceleration due to a cosmological constant or not? In many ways, an affirmative answer 
is the most depressing, since a true cosmological constant varies in neither space nor 
time, and thus, if it is the driver of acceleration, we already know all that we will ever 
know about it. As such, we would not expect further observational insights into the 
required unnaturally small value. More generally, we would like to be able to distinguish 
among this minimal explanation, a possible dark energy component and the option of 
an infrared modification of gravity. It has become increasingly clear that our best hope 
for this is to compare geometrical measures, such as the distance to the surface of last 
scattering, with measures of the growth of large-scale structure, which depend in detail, 
in a redshift-dependent manner on the gravitational force law underlying the collapse 
of overdensities. 

In this article we have endeavored to provide an overview of the major approaches 
to cosmic acceleration and to the techniques through which we hope to further our 
understanding. What should be clear is that the issues involved here are far from 
settled, and that the field is hungry for new ideas. It is unclear where such an idea will 
come from, and for this reason it seems important to pursue all possible avenues, no 
matter how unpromising they may seem. We hope to have provided a snapshot of some 
of these avenues of enquiry. It remains to be seen which, if any, of these will provide 
the key to the problem of cosmic acceleration. 
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